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Program No: FM - 110 
Orig. Date: 1 Dec 59 
Programmer: H. R. Knitter 

Program Title : Random Number Generator (2 methods). 

Purpose : To provide a basis for random choice within a program. 

Method 1 : A set of transformations are performed on a given number N so that 
the resultant N 1 is sufficiently "unpredictable", A double length 
number N is formed, then N* = (N '**") least sig. bits plus 
(N hS ) most sig. bits where overflow is ignored. 

Procedure : The following are the necessary orders: 

SQR n 

MPY n 

STO n 

XAR - 

ADD n 

TOV next instr. 

STO n 
Where n is the name of the location for N and where the 
resultant N 1 will be found. 

Output : The output is a nuriber 0<N'< 1 at binary 0. 

Limitations : This method has a reasonably large cycle time for all 
numbers tested as the original N. Good results were 
obtained using +10.^732. 0+2^, 0U6l,0 which had a cycle 
greater than ^5,000 repetitions. This method 
degenerates when N = but this fortunately, seems to 
generate in only rare cases. This method is not 
recommended for advanced scientific determinations, 
(e.g. using Monte Carlo methods) since the only test 
other than cycle length, has been distribution, which 
was sufficiently square between and 1. 

Note : To obtain a random number A such that 0^A<B where B is a 
given maximum, multiply the original random number N 
(at binary 0) by B. The resultant number A will be 
random and have the same scaling as B, 

Method 2 : The method used is the one suggested in the Journal of the Assoc, 
for Computing Mach. , Oct 59> page 527. The method essentially 
defines a new random number Nj as [N^., + N^l^^at binary 0. 

Procedure : A transfer to the initial location (i^ ) will initiate 
the generation of a new random number which will be 
in the accumulator upon return. 

Note: Same taote as in method 1 above. 



Because of method l f s short length and simplicity, a program tape for this 
method will not be distributed. The number FM-110 refers to method 2. 
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67HO.O +15.7760.0 .SAX 
+64.67*K).Ol'XTL 

67^1 *o +66.6750.0 Vctv 

+57.7(62.0 TRA 

67*12.0 +01 .7773.0 ADD V&7 

+42.7772.1 STA S+ir* 4er«r* 

67*13.0 +00.7777.0 cla /? - //3> /?**;//*" 
+01.7774.0 ADD /p f„ fat- 

67^.0+^2.7771.0 STA (Jim.) 
+53.6000.0 TOV 

67^5.0 +00.7777.0 CLA f 

+OI.7776.O ADD Jjt*,3r 

67*6.0 +33.7775*0 7 EXT a m .J 

+OO.6757.O VSTO „ * 

67^7.0+01.777^.0 ADD A, 

+1*2.7770.1 STA Ut) 

675O.O +H2.7772.0 STA Uk } 
+00.0000.0 CLA 1 

675I.O+OI.OOOO.O ADD \&&fi*J +[/bJ m 0jjL H 
+53.7772.0 TOV ) J u "'« (ft 

6752.O +60.0000.0 STO 

+57{7765{5f^TRA — £-x/r 

6753.0 +00.0000.0 . ,. 

-00.0000.1 V * 

6754.0 +00.6760.0 — 

-OO.6760.O — Ao 

6755.0 +00.0017.0 — , 

-00.0017.0 -p.- W*ff< 

6756.O +00.0001.0 t „-*<* 

-00.0001 .0 — /tf/T, 3r 

6757.0 +00.000*1.0 

-00.0004.0 — /; 
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6760.0 
6761 .0 
6762.0 
6763.0 
676^.0 
6765.0 
6766.0 
6767.0 

6770.0 
6771.0 
6772.0 

6773.0 

677^.0 

6775.0 

6776.0 
6777.0 



^•?137.o 

+06.2665.0 
+61.1755. 

+H3.3362. 
+16.3662. 

+11.675*1. 

+11.7061. 



+72.2703 # 
+IO.O655. 



+32.5326.0 



+55. 
+70I00I 

+07.1^01. 
-O6.6127. 



+01 .5200. 

-O5.7266. 
+11.1550* 

•^7*3735. 



.0 — 



+66.' 
+35-i 



a 



^365*0 

+66.W56.0 
-23.2217.0 

+1*0.761^.0 
-16.3650.1 

+50.6361 .0 
+65.3027.1 

+05.3301.0 
-2U73050O 
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